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1. Introduction 



The number of colours N c is an implicit parameter of QCD. In the large N c limit the theory 
becomes simpler amongst other features, mesons and glueballs become exactly stable. The 
finite-A^ corrections in the pure gauge theory are ff(l/N 2 ). This is also the case in quenched 
lattice field theory, since in that case fermion loops do not contribute. For full QCD the corrections 
are 0(1 /N c ). Clearly it would be very interesting to evaluate these corrections and establish how 
large they are in the case of QCD (N c = 3). 

The corrections can be estimated by carrying out quenched calculations at finite N c . The lead- 
ing N c dependence will then give the &(\jN^) quenched corrections, and an N c — > °° extrapolation 
will give the large N c limit. Finally this can be compared to unquenched SU (3) results (or indeed 
directly to experiment) to estimate the &(l/N c ) unquenched corrections. 

Calculation of quantities in the large-A^. limit is also important for attempts to extend the 
AdS/CFT correspondence [Q] to the case of QCD. This approach, known as AdS/QCD, attempts 
to learn about QCD by studying a five-dimensional theory that is dual to the large-A^. limit of 
QCD [Q, |], |5j, ^, ||]. The form of this five-dimensional dual (if it exists) is unknown, and 
additional information on the QCD side would help to constrain it. Since the duality occurs in the 
large-A^ limit such information must also be calculated in that limit. 

The masses of low-lying glueballs have been calculated on the lattice for N c up to 8 0. The 
finite-^ corrections were found to be small, reaching only about 10% in the case of SU (2). The 
pion mass has been calculated for N c = 17 to 23, with no N c dependence observed Jlo|]. Other 
meson masses have not to our knowledege been calculated at large N c . 

In this work we present calculations of the pion and rho masses up to SU (6). We describe our 
methods in the next Section and our results in Section [| We sum up in Section [| 



2. Methods 



We use the lattice software package Chroma [|1 1|], which we adapt to work for arbitrary N c . We 



use unimproved Wilson fermions, together with the Wilson plaquette action for the gauge fields. 



To set the scale, we use the string tension calculations by Lucini et al. []12|]. We choose the 
coupling j6 = 2N c /g 2 = 2Af?/A, where A is the 't Hooft coupling, such that the string tension in 
lattice units, a^fc, is the same for each A^. We use the value a\fc = 0.2093: for SU (3) this 
corresponds to j3 = 6.0175. The values for other Ac are shown in Table [I]. Adopting a value of 420 
MeV for the string tension, the lattice spacing is then 0.099 fm in each case. The values of a^fo 
used in the fits in []l2|] are very accurate so the errors on our estimate of the lattice spacing are 
less than 1%. 

The volume of our lattices is 16 3 x 32 in lattice units, corresponding to a spatial volume of 
1.58 3 fm. Finite-volume effects are expected to decrease with N c , and to be zero at infinite as 



long as the box is larger than a critical length l c [13]. This means that we should obtain the correct 



large-A^ limits for the masses we calculate, despite any finite-volume effects at small N c . 

We find that, as expected, the cost of updating a gauge configuration is approximately propor- 
tional to A" 3 , and the cost of inverting a propagator is approximately proportional to N% (the number 
of conjugate gradient steps required for the inversion is approximately independent of N c ). For the 
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N c 


P 


2 


2.4645 


3 


6.0175 


4 


11.028 


6 


25.452 



Table 1: Value of the coupling j3 for each N c . 

relatively small values of N c we use, the latter in fact dominates. Furthermore, we find that the 
correlators become less noisy at larger N c (see below), so fewer are needed to calculate a mass to a 
given precision. This reduces the total cost to ~ N c . A heuristic argument supporting this observed 
reduction in noise is the increase of the degrees of freedom of the statistical system <=c jV~ at fixed 
volume. 

We calculate local-local and local-smeared correlators on our configurations. We smear using 
Oto 100 iterations of Gaussian smearing, y/(x) = \jf(x) + K£„(£/ jU (jt)v't J < : +M) + U$,(x — ix)\jf(x— 
/i)), with smearing parameter fc = 4. We smear the gauge links using 10 iterations of APE smear- 
ing, U'p (x) = aU ti (x) + '£ v (U v (x)U^(x + v)ul(x + n) + ul(x-v)U fl (x-v)U v (x + n-v)), with 
smearing parameter a = 2.5 

3. Results 

We search for the critical value of the hopping parameter, Kc, by finding the value of k for 
which the pion mass m n vanishes. We show a plot (a«%) 2 as a function of 1/jc in Fig. [j] For 
suffciently small masses the dependence is linear, as expected, and we can extrapolate to zero pion 
mass. The values of K c we obtain are shown in Table ||. It is clear that they are rapidly converging 
to an 7Y C = °° value. 



N c 




2 


0.15327(16) 


3 


0.156397(45) 


4 


0.158168(39) 


6 


0.159060(44) 



Table 2: Critical hopping parameter K c for each N c - 

As mentioned above, we find that correlators become less noisy as A^ c increases. We illustrate 
this in Fig. where we compare pseudoscalar correlators on individual gauge configurations in 
St/ (3) and SU(6). The pion masses are almost identical, am n = 0.283(5) and 0.286(3) respec- 
tively, but the scatter between individual configurations, a measure of the noise, is about twice as 
large in SU (3), in agreement with the naive degrees of freedom argument above. 

Since fluctuations decrease as A^ c increases, and we use correlators of fluctations to extract 
masses, one might worry that we will be unable to calculate masses in this way in the large-Af c limit. 
Fortunately this is not the case. The correlators do indeed decrease as 1/iVf, but the fluctations in 
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Figure 1: (a»%) 2 as a function of 1/fC for 5C/(2) (+), SU(3) (x), 5?7(4) (*) and St/ (6) (□). 
10 




Figure 2: Point-point pseudoscalar correlators on individual gauge configurations in 51/(3) at K" 
0.1547 (+) and in SU (6) at K = 0.15715 (x). SU(6) results have been shifted vertically for clarity. 



the correlators decrease at the same rate, so the signal-to-noise ratio remains constant. For a more 
detailed discussion of this issue see 

We have calculated the lowest Dirac eigenvalues at our lightest pion mass. We plot the distri- 
butions of eigenvalues we obtain in SU (3) and SU (6) in Fig. ||. These are at similar pion masses, 
am % = 0.188(6) and 0.212(3) respectively, but we see that the distribution of eigenvalues is much 
narrower in SU(6), with no near-exceptional configurations. This suggests that it may be possi- 
ble to use unimproved Wilson quarks at much lighter pion masses at large N c than is possible for 
SU{3). 
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Figure 3: Eigenvalue density for the smallest eigenvalue of 75D in SU (3) at K — 0.1557 (red) and in SU (6) 
at k: = 0.1581 (green). 



We extract the ground state pseudoscalar and vector masses from correlators at a range of 
quark masses. The lowest pion mass is approximately m K = O.laT 1 =410 MeV, corresponding 
to a quark mass somewhat lighter than the strange quark mass, while the heaviest is around m K = 
0.6a _1 = 1220 MeV. At our lighest masses m n jm, p « 0.5. We show the extracted masses in Fig. |j. 
We do not have exact matching of the pion masses for different N c in all cases, so we plot the rho 
masses against the squares of the pion masses — if there is no 7V c -dependence these should all fall 
on a common curve. We see that for the entire range of quark masses, the points do indeed fall 
onto the same line, the only exception being SU (2), for which m p appears to be about 10% higher 
at the lowest quark masses. Finite volume effects are expected to be largest at small N c and light 
quark masses, but we would expect these to increase m\ by a larger amount than m p , shifting the 
SU (2) points below the line, the wrong direction to explain the observations. If the deviations are 
indeed real finite-Af c effects, they are of order 10%, similar to the deviations observed for glueballs 
in Assuming they are dominated by the l/Ng corrections, the corresponding deviations for 
St/ (3) should be around 4%. 



4. Conclusions 

We have calculated the masses of the pion and rho mesons up to SU (6), at a cost approximately 
proportional to N c . We have shown that the N c dependence of the masses is small, approximately 
10% for SU (2) at small pion masses, and that there is no observable dependence at higher masses. 
If these differences are dominated by the 1 /N 2 corrections, which seems very plausible, there will 
be a difference of about 4% between SU (3) and SU (°°). Thus, at least for the masses of these 
mesons, quenched QCD is very close to the large-Af c limit. 

The distribution of the lowest eigenvalue of 75D becomes much narrower as N c increases. This 
suggests that it maybe be possible to use Wilson quarks at significantly lower quark masses without 
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Figure 4: Rho mass against pion mass squared, for St/ (2) (+), SU(3) (x), SU(4) (*) and St/ (6) (□). 



running into exceptional configurations. We intend to explore this possibility in future work. We 
will also examine smaller volumes to check that the volume dependence does indeed decrease. 

Apart from the ground state pseudoscalar and vector, we intend to also calculate the masses of 
the excited states. It would also be interesting to look at other quantum numbers, in particular the 
scalar mesons. 

Finally we note that our calculations have been carried out at only one lattice spacing. Our re- 
sults will thus be affected by lattice corrections. However, these lattice corrections will themselves 
have a large-Af c limit, which will affect our results for all N c equally. Only an unlikely cancel- 
lation between the continuum and lattice 1 /N% corrections could cause us to see the very small 
/^-dependence we observe. Thus, although our values for the meson masses have yet to be extrap- 
olated to the continuum limit, our conclusion that the -dependence is very small is unlikely to be 
affected. One should of course check this by repeating the calculations at a smaller lattice spacing, 
and we intend to do so in future. 
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